arXiv:gr-qc/9507043vl 20 Jul 1995 


A Model for Stars of Interacting Bosons and Fermions 


Claudio M. G. de Sousa [| and J. L. Tomazelli | 

International Centre for Thoretical Physics, P.O. Box 586, 34100 

Trieste, Italy 


Abstract 


In this paper we introduce a current-current type interaction term in the Lagrangian density 
of gravity coupled to complex scalar fields, in the presence of a degenerated Fermi gas. For low 
transferred momenta such a term, which might account for the interaction among boson and fermion 
constituents of compact stellar objects, is subsequently reduced to a quadratic one in the scalar 
sector. This procedure enforces the use of a complex radial field counterpart in the equations of 
motion. The real and the imaginary components of the scalar field exhibit different behaviour as 
the interaction increases. The results also suggest that the Bose-Fermi system undergoes a BCS-like 
phase transition for a suitable choice of the coupling constant. 
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1. Introduction 


In modern approaches to astrophysics and cosmology, concepts that were initialy re¬ 
stricted to elementary particle physics acquire a wider meaning, particularly in inflationary 
models. These models give predictions for the mass density of the present Universe larger 
than that observed, if one presumes that it is close to the critical value. This suggests that 
there is a large amount of hidden matter, which has not been detected so far. 

Among the possible candidates for the so called dark matter are boson stars^^l, which 
consist of gravitational bound states of scalar particles. Possibly, such structures were formed 
through gravitational collapse in the early Universe and may appear also in the core of 
composite objects, whose external envelope are made of standard baryonic matter. 

Since in addition to the bosons there were also fermions in the primordial gas, we would 
expect boson-fermion stars to prevail. This system was studied in detail by Henriques et al 
[2]. In their work it has been shown that the properties of boson-fermion stars are qualita¬ 
tively the same, irrespective of the addition of a self-coupling term for bosons. Nevertheless, 
it seems that there is no modeling for explicitly dealing with the interaction between bosons 
and fermions in such systems. 

In this work we introduce an effective coupling between bosons and fermions to afford 
a more realistic description of the system and compare our results with the ones in the 
current literature. In section 2 we construct the energy-momentum tensor for interacting 
fermions and bosons using the Schwarzschild metric. In section 3 we obtain the evolution 
equations for the coefficients of the metric and for the fields. In section 4 we exhibit the 
results of numerical simulations for the corresponding dynamical system. Finally, in section 
5 we discuss our results and compare with those obtained without taking the interaction 
into account. 


2. The Boson-fermion Interaction 


Before introducing the interaction term, we outline the non-interactive boson-fermion 
model. We assume the metric to be the standard Schwarzschild one, 

= —B{r)dt^ + A{r)dr‘^ + r'^dd'^ + sin^ 6d(p‘^ . (1) 

Our sign and conventions are the same as those used by Fiddle The Lagrangian for 
the complex scalar field with no explicit interaction is given by 

C = , (2) 

IGttG 

where 

<I>(r, r) = (3) 

and for our purposes 0 is a complex scalar field. 

For the combined bosons-fermions star we consider that the fermions are described by a 
perfect fluid with energy density p and pressure p as proposed by Chandrasekhar dl; 
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p = i^(sinht — t) , 

K t 

p = —(sinht —8sinh—h 3t) , 
3 2 


( 4 ) 


where t is a parameter, K = mf^/32n^, and is the fermion mass (to be considered as that 
of neutrons for illustrative purposes). The evolution equation for the fermions is given by 
an equation of state, namely we have 


P = -^(P + P)^ ’ 


( 5 ) 


where the primes stand for derivatives whith respect to r. 

The corresponding energy-momentum tensor for the bosons and fermions without inter¬ 
action reads as 


7^(0) _ rpB I rpF 

flip flip ' flip 


( 6 ) 


with 


= (P + P)u^u^ + P9t,u , (8) 

where superscripts B and F label bosons and fermions from now on. 

At this point we introduce the following interaction term in the Lagrangian density 

, (9) 

where 

J^($) = , (10) 

f(pP) = (11) 


which represent the boson and fermion currents, respectively, while the y’s are the usual 
Dirac matrices, which satisfy 7°”^ = —7° and 7*^ = 7® (i = 1,2,3). This is a typical contact 
or current-current interaction between bosons and fermions where the coupling constant has 
dimension [A] = M“^. We emphasize that in dealing with a nonrenormalizable interaction 
term an energy scale must be introduced when we consider quantum corrections; we hnd a 
similar situation in the pure Einstein-Hilbert gravity. Note that this point is of no relevance 
since we are giving a semiclassical treatment to the problem. 

The complete Lagrangian density, including the interaction term, is invariant under global 
U{1) gauge transformations. For low transferred momenta we may replace Eq.(8) by 

fW = , (12) 

where 

r° = iu° , r = u^. (13) 
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The four-vector = (m°, Ur, ug, u^) is the four-velocity of the fermion fluid. 

The contribution of the interaction term for the energy-momentum tensor is given by 

T;f = , (14) 

which together with (6) give the total energy-momentum tensor 

rji _ ^0 I rpint 

~ fiiy ' fiu * 

Using (3) and its complex conjugate in (14) we arrive at 


T»“ = -■pr'm, 

Tp”' = -2u^4>‘<»i’■ 


(15) 

(16) 


Instead of the bilinear operator ipip we take the ground-state of the fermions system 
< ipip >F in the semi-classical approximation. In this limit we replace this quantity by the 
average fermion density np. 

Next, we write the complex scalar helds as 


0 — 01 + *02 , 0 * — 01 — *02 , 

which yields 

^ ^ (17) 

= -^(0? + 0^), (18) 

where a = 2nFX. The remaining nonvanishing contributions for the total energy momentum- 
tensor follows immediately from (6). 


3. Evolution Equations 


Now we are ready to write the evolution equations for our system. With an appropriate 
redehnition of the dynamical variables and parameters, namely 


X = 

mr , 

a{x) = 

V87rG0(r 

47rG 

W) = 

AtiG 

pit) = 


a 

u 

a = — 
m 

, w = — 
m 
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the equations read 


A' 


B' 


0-1 

(^2 


Si 


S2 


f 


t' 


xA^ 

xAB 

Si , 

, 


W 


a 




X 


f \ 9 s"^ wa ,1 ^ , 

2p+ _-q^ + --_.VT(-4-l). 


B 

X 


A{^-l + ^]a, + ai/As,+ 


aw 


B 


Vb 


a 


^/A ( B' 2' 




2B ' XI ^ ' 


-A 


w 


aw 


B ^/B, 


0-2 


— ay/Asi + 


1 /A' _ B'" 
2\'A~^ 


1 fA' _ B'' 

2 vl" ~ ^ 


Si 


S2 


-2' 


ay/A ( B' 2\ 

2 ^ xj ’ 

R/ sinh t — 2 sinh - 

^_ 2 _ 

R t 

cosh t — 4 cosh 2+2 


(19) 

( 20 ) 

( 21 ) 

( 22 ) 


(23) 


(24) 

(25) 


where 

2 2,2 2 2,2 

(j — (jY A (^2 1 s — -|- S 2 • 

These equations form a non-autonomous system of non-linear hrst-order differential equa¬ 
tions, which cannot be linearized due to the quadratic terms involved. Equations (19) and 
(20) are the Einstein equations; equations (21)-(24) correspond to the Klein-Gordon equa¬ 
tion, whereas (25) gives the evolution of the fermion energy density and pressure. In the 
above set of equations the dynamical variables and parameters are dimensionless. From now 
on the primes stand for derivatives with respect to x. 

It is noteworthy that these equations are invariant under the scale transformation 


B —rjB , w — y/^w , 

even after the inclusion of the interaction term. Since the initial value Bq is undetermined, 
this permits its redehnition during numerical calculations, in such a way to obtain the 
asymptotic values for the metric coefficients converging to those of a flat space. 


4. Numerical Results 


In this section we present the numerical simulation results of the above set of equations. 
The set of equations were solved by using the hfth-order Runge-Kutta method; we also use 
an apropriate “shooting method” to infer the value of w, according to an initial value of Bq. 
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For further details on the numerical criteria, the interested reader is referred to Rufhni and 
Bonazzola^^l We require that the metric to be asymptotically flat and that the scalar helds 
as well as their derivatives to vanish at inhnity. 

Close to the singularity at the origin, which is inherent to Schwarzschild metric, we see 
that the last terms in (19) and (20) are dominant. In this region A ^ xjix — const.) and 
B ^ [x — const.)/x, so that B and B' diverges and A becomes oscillating. Hence, for 
numerical purposes it is convenient to start with x slightly shifted from the origin. 

In Fig.l we show the results hrst obtained in Ref. [5] for a typical boson star, where we 
note that there is no overlap between a and a', which characterizes the ground state of the 
system. Throughout the simulations the initial values are: = 0.32, s? = = 0 and 

Aq = 1. 

In Fig.2 we show the ground-state of a boson-fermion star for to = 7. In this case 
Bq = 0.036 and w = 1.098. 

Fig.3 displays the same curves for to = 8.7 which corresponds to Bq = 0.01 and w = 1.327. 
Notice that the peak in the curve for A has increased and has been shifted towards the 
origin, while a and a' approach to zero faster. Larger values of to mean higher fermion 
energy densities, so that the fermion contribution to the energy-momentum tensor becomes 
dominant. In this sense we would expect curves like those given in Fig.4 to be in agreement 
with the pioneering results of Oppenheimer and Volkoff for neutron stars. 

Figures 5 and 6 exhibit the results after introducting the interaction term, for the choice 
to = 8.7, and different values of the dimensionless coupling constant a. When we switch the 
interaction on at small values of a, w increases, as shown in table 1. On the other hand t 
vanishes at a smaller x. 

If we continue to increase the interaction we observe that, for a certain value of a, w 
suddenly decreases, suggesting that there is a critical value of a, in the range 10“^ — 10“^, 
in which the system experiences a second-order phase transition, corresponding to boson- 
fermion pair formation. This phenomenum is similar to that observed in the BCS theory; 
this is not so surprising since we have started with a quartic interaction term of the form 
(9). In Fig.6 (b) and (c) we can observe the splitting of the real and imaginary parts of the 
scalar held that also occurs in such interval. For completeness, we also exhibit the phase 
space diagrams for the metric coefficients and for the scalar helds at a = 10“^ in Fig. 7. 


a 

w 

Bo 

10-1^ 

1.482 

0.11 

10“^ 

1.557 

0.11 

10“® 

1.565 

0.011 

10“=^ 

1.335 

0.010 

10-2 

1.060 

0.009 


Table 1: values of w an Bo for different coupling constants a and to = 8.7. 
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5. Concluding Remarks 


In this work we studied a model for boson-fermion stars in interaction, in the region 
of low frequencies, and observed the behaviour of the system for increasing values of the 
coupling constant a. At small values of a there is no signihcant changes in the system; as a 
becomes larger and larger the ground state energy of the system increases and the fermion 
energy density and pressure vanish at smaller distances. As a result, the fermion star is 
conhned to a smaller region, after which the scalar helds are still present. 

When a reaches a critical value, the ground state energy of the system suddenly decreases, 
indicating the possible occurence of a second-order phase transition. In this case, we might 
have supplied enough energy to bind bosons and fermions together. 

From the results outlined above it is possible to compute the mass and the radius for 
typical boson-fermions stars; this is the subject of a forthcoming paper. Another interesting 
matter, is to develop a qualitative approach to our evolution equations with Schwarzschild 
metric, in order to compare the phase diagrams with our corresponding numerical results. 

A more general treatment of boson-fermion gravitationally bounded systems should incor¬ 
porate spin effects by considering the full Dirac Lagrangian. However, it would be necessary 
to extend the system of differential equations from seven to hfteen, which might demand a 
great effort. 
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Figures Captions 


Figure 1: The values of the helds and metric coefficients in a typical boson star conhgura- 
tion. The corresponding value of w is around 0.826. 

Figure 2: Fields and metric coefficients for a typical boson-fermion star, with R = 7. In 
this case, the value of w has increased to 1.098. 

Figure 3: The same curves as those of Fig.2, for fg = 8-7. 

Figure 4: The evolution of t for different initial values. 

Figure 5: Fields and metric coefficients for a boson-fermion star when we switch on the 
interaction at a = 10“®. 

Figure 6: The same curves as those shown in the last hgure, for a = 10“^. 

Figure 7: Phase space diagrams (a) for the metric coefficients and (b) for the scalar fields. 



